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We study the strangeness contribution to nucleon matrix elements using Nf - = 2 + 1 dynamical 
clover fermion configurations generated by the CP-PACS/JLQCD collaboration. In order to eval- 
uate the disconnected insertion (DI), we use the Z(4) stochastic method, along with unbiased 
subtraction from the hopping parameter expansion which reduces the off-diagonal noises in the 
stochastic method. Furthermore, we find that using many nucleon sources for each configuration 
is effective in improving the signal. Our results for the quark contribution to the first moment 
(x) q in the DI, and the strangeness magnetic moment show that the statistical errors are under 
control with these techniques. We also study the gluonic contribution to the nucleon using the 
overlap operator to construct the gauge field tensor, F^ v . The application to the calculation of first 
moment, (x)g, gives a good signal in quenched lattice QCD. 
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1. Introduction 

The understanding of the structure of nucleon is one of the central issues in hadron physics. 
For instance, the parton distribution functions (PDFs) have been studied extensively, and the ob- 
servation of scaling violation in PDFs provides the cradle for the establishment of the fundamental 
theory, QCD. Yet, there exist many unresolved questions for the structure of the nucleon. The EMC 
experiments [|l|] and subsequent experiments show that quark spin carries only ~ 30% of the total 
spin of the nucleon. Consequently, one concludes that the remaining ~ 70% should be carried by 
the quark orbital contribution and the glue. Direct quantitative identification of each contribution 
has been undertaken by experiments. The strangeness contribution to the nucleon structure is also 
under intensive study experimentally. 

Under these circumstances, it is desirable to provide definitive quantitative results using the 
lattice QCD method. In fact, there are many lattice QCD studies of the nucleon structure [Q]. 
However, most of these calculations are limited to the so-called "connected insertion (CI)", and 
there have been few calculations considering "disconnected insertion (DI)" [||, |], ||]. Although 
the calculation of DI is known to be a very difficult problem, we emphasize that DI are related to 
rich physics, e.g., only DI consists of the the strangeness contribution to the nucleon. We shall 
describe our methodology to obtain the signal effectively in DI calculation. We also note that using 
Nf = 2 + 1 dynamical configurations could be essential for small quark masses. In this proceeding, 
we present the study of the DI part for the quark contribution to the first moment of the nucleon, 
{x) q , whereas the study of the CI part is presented elsewhere [§]. 

Another very important, but often not calculated component, is the glue contribution to the 
nucleon structure. This is because the straight-forward calculation using the standard link variables 
is known to yield very noisy signal We have proposed [[|] to use the overlap operator (D ov ) 
to overcome this problem, and we show the first application of Tr s [a^ v D ov ] to the gluonic first 
moment, (x}g- 

2. Formalism and simulation parameters 

We use the Nf = 2 + 1 dynamical clover fermion with renormalization group improved gauge 
configurations generated by CP-PACS/JLQCD collaboration [g]. We use the j3 = 1.83 configura- 
tions with the lattice size of V = L 3 x T = 16 3 x 32, for which the lattice unit is a -1 = 1.62GeV 
and physical spatial size is (2fm) 3 . The hopping parameters for light (u,d) quarks are K = 0.13825, 
0.13800 and 0.13760, which correspond to m n = 610, 700 and 840 MeV, respectively, and the 
hopping parameter for strange quark is fixed to be K s = 0.13760. We perform the calculation only 
for the dynamical quark mass points. For each quark masses, about 800 configurations are used. 

We also perform the complementary study using the Wilson fermion with Wilson gauge action 
in the quenched approximation. We generate 500 configurations of 16 x 24 lattice at /3 = 6.0, 
where the lattice spacing is aT x = 1.74GeV and the physical spatial size is (1.8fm) 3 [Eh. The 
calculation is performed with three light quark hopping parameters of K = 0.1540,0.1550,0.1555, 
which correspond to m % = 480 — 650MeV, with strange quark hopping parameter fixed to K s = 
0.1540. 
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The nucleon matrix elements can be calculated by taking the ratio of three point function 
to two point function n 2pt , 

Uf(p,t 2 ; q,h; P,x ,t ) = Y,J~ ip '* 2e+q '*HQ\T[JN(Z2j2)0(xiA)JN{xo,to)}\O), (2.1) 

X 2 ,Xl 

U 2 v\p,t; x ,t ) = £e-V*(0\T[Jiffct)JN&,t )] |0>, (2.2) 

where G is an appropriate operator for the matrix element of concern and is the nucleon inter- 
polating field. 

The calculations of three point functions of DI involve evaluation of both of the two point 
function part and the quark loop part. Because the latter requires all-to-all propagators, for which 
the straightforward calculations are practically impossible, we use the stochastic method [[K]] as 
follows 




(2-3) 



where F is an arbitrary matrix and tj corresponds to the l-th noise. In the practical calculation, 
this method introduces a variance, because L = N no i se is finite. In order to reduce such off-diagonal 
error [JHJ, we use the unbiased subtraction from the hopping parameter expansion (HPE) [p"T|]. 
In our practical calculation, we adopt Z(4) noises in color, spin and space-time indices. We take 
Nnoise = 300(500) for each configuration for full (quenched) QCD simulation, respectively, along 
with the use of HPE up to the 4th order. 

In the stochastic method, it is quite expensive to achieve a good signal to noise ratio (S/N) by 
just increasing N, w j se because S/N improves with \/N no i se . In view of this, we use many nucleon 
sources N src in the evaluation of the two point function part for each configuration. Because the 
calculations of quark loop and those of two point functions are independent, this is expected to be 
an efficient way to increase statistics. In fact, as we will show later, we observe that S/N improves 
almost ideally, i.e., by a factor of \/N src . 

For the study of glue contribution to the nucleon structure, it is essential to find a suitable glue 
operator. In fact, glue operators constructed from link variables suffer from large fluctuations in 
high-frequency modes, which causes poor S/N in the calculation [^]. We propose [||] to use the 
gauge field tensor constructed from the overlap operator D ov as 

Ffiv(x) = const, x Tr s [cfy v D 0V (x,x)], (2.4) 

where Tr s corresponds to the trace in spinor space. The advantage of this formulation is that the 
ultraviolet fluctuation is expected to be suppressed due to the exponential local nature of D ov . In 
order to estimate D m (x,x), we again use the stochastic method. In this case, we treat the color 
and spin indices exactly and space-time indices are diluted for two sites separation on top of the 
even/odd dilution. Therefore, the minimal length to the next neighbor site amounts to four hoppings 
away ("taxi driver distance"=4). We use two Z(4) noises and take the average between them for 
each configuration. 
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3. First moment of parton distribution 

Quark contribution to the first moment of the parton distribution in the nucleon, (x) q , can be 



obtained by using the following energy-momentum tensor operator [|1 2fJ , 
and by taking the following ratio of three point to two point function, 



T 4i = (-1)*- 



(3.1) 



Tr 



r e-^{p,h; (Mi ; p,xoM 



x e 



-iqx 



(2 Pi ) x Tr[r e •n 2 P t (p,f 2 ; *b,'o)] x <? +/? 



■V() 



(3.2) 



where r, 



I 





in spinor space (for the Dirac representation.) In order to improve the S/N, 



we further take the summation for the operator insertion time t\ for the range t\ = [to + — 1], 
where ?o( f 2) is the nucleon source (sink) time, respectively. 

In Fig. p], we plot the ratio of three point to two point functions for the strangeness, (x) s , in 
terms of t%. Because of the summation of operator insertion time t\, the linear slope corresponds to 
(x) s . Blue points denote the result for N src = 1 and red points for N src = 32. One can clearly see 
that increasing N src reduces the error significantly (by about a factor of \/N src ), and a clear signal 
can be extracted from the N src = 32 data. 

In Fig. ||, we plot the ratio of (x) s / (x) ut j(DI) for each quark mass. By taking the linear chiral 
extrapolation in terms of m 2 n , we obtain a preliminary result 



(x) s /(x) Utd (DI) =0.857(40), 



(3.3) 



in the chiral limit. Although it is necessary to consider the renormalization factor, this result shows 
that the statistical error is well under control and reliable lattice QCD calculation is possible even 
for the DI. 
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Figure 1: The ratio of three point to two point func- Figure 2: Ratio (x) s /{x) ud (DI) plotted in terms of 

tions plotted in terms of the nucleon sink time t 2 for m \ for dynamical configurations. Blue line corre- 

dynamical configurations with K = 0.13760. The sponds to the linear chiral extrapolation, 
linear slope corresponds to (x) s . 
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We perform the same calculation for the quenched case, and the result using N src = 16 with 



perturbative renormalization is Q13Q 



{x) s /{x) u , d (DI) = 0.88(7). 



(3-4) 



Gluonic contribution to the first moment, (x)g, can be obtained in the same way, by using the 
following energy-momentum tensor 



T 4 k = - £(+0 * (F 4i F ki ), 
i=i 



(3.5) 



where F^ v is constructed using the overlap operator as shown in Eq. ( |2.4[ ). In Fig. g, we plot the 
ratio of three point to two point functions in terms of the nucleon sink time ?2, using the quenched 
QCD configurations. As can be seen, we obtain a prominent linear signal for (x)g- It is quite 
encouraging that we obtain the signal with about three sigma accuracy. With unbiased subtraction, 
this could be improved to larger than four sigma accuracy. We are also working on the deflation 
which could improve the S/N substantially. Further investigation including the renormalization 
factor is in progress. 



=0.15400{utl), 0.15400(s) 



w/o subt (#src-l(i) ■ 



1 -Slope before subtraction = 0.069(26) 



, i I l l 
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I Summation of Insertion Time: (source time+1) to (sink time-1) 



Figure 3: The ratio of three point to two point functions plotted in terms of nucleon sink time ?2 for quenched 
configurations with K = 0.1540. The linear slope corresponds to (x)g- 



4. Strangeness electric and magnetic form factors 

We use the point-split vector current as 

= q(x) ( 1 - 7 M )Up {x)q{x + jU ) - q(x + JU ) ( 1 + Yfi ) C/J (x)q(x) . 



(4.1) 



The advantage of the point-split operator is that there is no additional renormalization factor to be 
considered, because it is a conserved current. Using this operator, the electric and magnetic form 
factors can be obtained ^ by the following combination of three point and two point functions: 



Tr 



xr» Tr 



r e -n 2 v%h; jfo,fo) 



Tr[r e -n 2 Pt(^?i; 4d))]x^ 



Tr 



r e -n2pt(S,* 2 ; x ,t ) 



E a — m 
2 ^ i F 



2m 



(4.2) 
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Tr 



r k -nX = .(o,t 2 ; qfy; -q,x Q ,t ) 



xe ->q*o Xr 



r e -n?i*(6,ti; xo,t ) 



Tr [r e • T&*(q, h ; x ,t )} x e+&*> 



Tr 



(-1): 



1 



E a + m 



e !7 ^(F 1 (-4)+F 2 (-^)), 



(4.3) 



where I\ = 



a* 



, m denotes nucleon mass, and E q = \/m 2 + q 2 . 



o o y ' * 

In Fig. |4|, we plot the ratio of three point to two point functions in terms of the nucleon sink 
time ?2 where the linear slope corresponds to the signal of strangeness magnetic form factor at 
q 2 = (2n/L) 2 . In order to obtain the final quantitative result for the strangeness magnetic moment, 
detailed study for several Q 2 (—q 2 ) and chiral extrapolation are needed and they are in progress. 



K^O.l^filKud), 0.137601s) 
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Figure 4: The ratio of three point to two point functions plotted in terms of the nucleon sink time tj_ for 
dynamical configurations with K = 0.13760. The linear slope corresponds to the strangeness magnetic form 
factor at q 2 = {2n/L) 2 



5. Summary 

We have studied the strangeness and gluonic contribution to the nucleon matrix elements from 
lattice QCD. The strangeness matrix elements have been studied by the Z(4) stochastic method, 
with the unbiased subtraction from the hopping parameter expansion in order to reduce the off- 
diagonal noises. We have further taken many nucleon sources for each configuration, and observed 
that this method is almost ideally effective to improve the signal with modest cost. Using Nf = 
2 + 1 dynamical clover fermion configurations, we have analyzed the quark contribution to the first 
moment, (x) q , and the strangeness magnetic form factor, and found that statistical errors are well 
under control with these improvements. 

The gluonic contribution for the first moment of nucleon, (x)g, has also been studied with 
the use of the overlap operator to construct the gauge field tensor, F^ v , and we have shown the 
effectiveness of this method by the explicit calculation at the quenched level. 

Although there remain several sources of the systematic error, such as the finite volume arti- 
fact, discretization error, excited-states contamination and chiral extrapolation, we plan to inves- 
tigate these issues by using the configurations with larger volume and lighter quark masses with 
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various lattice cut-offs. The study for other matrix elements such as the angular momentum contri- 
bution to the nucleon spin is also in progress. 
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